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Abstract 

The eigenvalue problem for the p-wave bound states formed by two unequal-mass 
scalar particles through the massive scalar particle exchange is analyzed numerically in 
the framework of the Bethe-Salpeter ladder model. As in the s-wave case, the eigenval- 
ues of the coupling constant are found to become complex for some mass configurations 
in some range of the bound state mass. The Bethe-Salpeter amplitudes of the low-lying 
bound states are also investigated. 
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1. Introduction 



The Bethe-Salpeter (BS) formalism is one of the effective (although now a rather old- 
fashioned) method to analyze the nonperturbative aspects in quantum field theory. ^ Usually, 
in the two-body homogeneous BS equation (describing two-particle bound states), complete 
propagators are approximated by free propagators and the BS kernel (describing the two- 
body interaction) is approximated by its Born term. For the case that the two-body inter- 
action is induced by a particle exchange, the approximated BS formalism is referred to as 
the ladder model. In this model, given a mass of the bound state, the BS equation can be 
interpreted as an eigenvalue equation for the coupling constant. 

Despite of these simplifying approximations, the two-body homogeneous BS equations 
cannot be treated analytically except for the Wick-Cutkosky model. ^ This model, describ- 
ing two scalar particles (with not neccessarily equal masses) exchanging a massless scalar 
particle, can be reduced to an eigenvalue equation of Sturm-Liouville type. It's spectral 
property is fully clarified analytically as well as numerically. In the case that the mass of 
the exchanged particle is nonzero (termed as scalar-scalar ladder model by Nakanishi 3 )), the 
eigenvalue problem remains, even after the partial wave decomposition, a two-dimensional 
partial differential (or integral) equation. If the two external scalar particles have equal 
masses, the corresponding eigenvalue equation can be reduced to a manifestly Hermitian 
form through the Wick rotation. In this case, although the detailed properties of eigenvalues 
are less well understood compared with the Wick-Cutkosky case, general properties like the 
discreteness, reality and positivity can be derived from the analysis of the integral kernel. 

Contrary to the equal-mass case above, in the unequal-mass scalar-scalar ladder model, 
even the reality of eigenvalues is not yet proved or disproved analytically. The situation 
was quite confusing nearly thirty years ago. There were coexisting the papers suggesting 
the reality of eigenvalues analytically 4 ^ and numerically 5 ^, and also the papers supporting 
the nonreality of eigenvalues analytically 6 ^ and numerically 7 ^. Relying on the progress in 
numerical computational technique, the present authors performed a rather extensive cal- 
culation 8 ) to obtain eigenvalues numerically. It was found that complex eigenvalues could 
appear in the first and second s-wave excited states, depending on the mass configuration 
of the two external particles and the exchanged particle, for the bound state mass squared 
around the pseudothreshold. The authors also made a numerical analysis of the BS ampli- 
tudes (eigenfunctions) 9 ) to investigate the behavior of amplitudes in the momentum space 
and also to check the numerical accuracy of the results obtained in Ref. 8. 

The purpose of the present paper is to study the "universality" of the complex eigen- 
values in the unequal-mass scalar-scalar ladder model based on the improved numerical 
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method compared with that of Ref. 8. As was pointed out in Ref. 8, a preliminary cal- 
culation suggested that all the eigenvalues of the p-wave bound states were real, for the 
mass configurations where certain eigenvalues of the s-wave states became complex. In the 
next section, the general formalism for the BS equation is presented. In §3, the eigenvalue 
equation in the p-wave case is solved numerically. It is found that if the mass of exchanged 
particle is lowered, complex eigenvalues appear also in the p-wave case. In §4, the radial 
parts of BS amplitudes of some typical bound states associated with complex as well as 
real eigenvalues are depicted in the momentum space. Physical meaning of the results is 
discussed and a further outlook is made in the final section. 



§2. General formalism 

The BS amplitude </>(p,P4) in the momentum space, with mass squared s (> 0), formed 
by two scalar particles with masses mi and m 2 by exchanging a scalar particle with mass /i 
in the ladder model, obeys, after the Wick rotation, the following BS equation: 

p 2 - (mVs + m) 2 ] [ml + p 2 - imVs-m) 2 <P(p,P4) 

A /"jo Hp 



~n 2 J P u 2 + (p-v') 2 + (v,-v' A ) 2 ' 1 ; 



/i 2 + (p-p') 2 + (p 4 ~p' A ) 

with r/i = mi/ {mi + m 2 ) and i]2 = m 2 /(mi + m 2 ). We shall denote the model as the 
[mi -<=/i=^m 2 ] model. The original setting of the bound state problem was to find the bound 
state mass y/s at a given value of the coupling strength A. Instead, due to the peculiarity of 
the ladder model, we can regard (2-1) as an eigenvalue equation for the coupling strength A 
at a given y/s in the range < \fs < mi + m 2 . Taking the unit mi + m 2 = 2 (A has mass 
dimension +2), we can transform Eq. (2-1), for the [(1 + A)<=/jl=^(1 — A)] model, into 

^ P,P4) = l(l + A) 2 (l-a)+p 2 -2t(l + A)^p A } 

x I A r 

[(l-A) 2 {l- ( r)+p 2 +2i(l-A)^p 4 \ 7i 2 J P ^ 2 + (p - p') 2 ' 1 j 

with a — s/4 and < A < 1. In this unit, the two-particle threshold is equal to 4, while the 

pseudothreshold is A A 2 . 

The three-dimensional (space part) rotational invariance of Eq. (2-2) assures that the 

four-variable integral equation can be reduced to a two-variable (e.g. p 4 and p := \Jp^ + p 2 ) 

equation. By using the polar coordinates in the four-dimensional momentum space (p, p^) , 

Pi = p cos (3, p = p sin /3(sin 9 cos ip, sin 9 sin if, cos 9) , 

(p > 0, 0<(3, 9 < 7T, < up < 2tt) 
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through the expansion of the BS amplitude in terms of four-dimensional spherical harmonics 
(smf3) l C l ] ^} l (cos{3)Yi m (9,(p) with associated normalization constants Nl,i, 

oo 

<t>(p,Pi) = X>L,i(p) • ^ • (sin/3)'Ci + _ 1 ,(cos/3)y im (^^) 

for suitable I (— 0, 1, 2, . . .) and m (— / <m<l), the two- variable equation is further trans- 
formed into an infinite system of one-variable integral equations. 

By changing the variable p (the magnitude of the four-dimensional momentum) to z 
defined by 

p=^—, z = ——^ (*:-l-l asp:0^oo) 
Eq. (2-2) is rewritten, after the partial wave decomposition, as 

oo r i 

g w ( z ) = Xj2 dz> •K LL ,, l (z,z')g L ,, l (z'). (L = l,l + 1,...) (2-3) 

L'=l J ~ X 

The Z-th wave BS amplitude 4>(p,Pa) is expressed in terms of Ql,i( z ) ' s as 

Mp,p,) = W 'g(ip) i - , ( T ^) w oK(^) »«g=i) n.(«,*>). (24) 

The detailed procedure of the above transformations and the explicit form of the integral 
kernel matrix functions K LL /j(z,z') for I = and I = 1 are given in Ref. 8. They are all 
expressed as finite sums of elementary functions of z and z' . These matrix elements are all 
real, but not symmetric (K L i Lj i(z' : z) ^ K LL ij(z,z')) in general. If Eq. (2-3) admit a real 
eigenvalue A, the corresponding eigenvector gL,i( z ) (L = 1,1 + 1, . . .) can be taken to be 
real. Even in this case, however, since imaginary unit i appears on the right-hand side of 
Eq. (2-4), the BS amplitude 4>{p,Pa) generally takes a complex value. 



§3. Eigenvalues 

In Ref. 8, we have analyzed numerically the eigenvalue problem (2-3) for the s-wave (I = 0) 
case. We have found there that, for the [1.6 <= 1.0 =^ 0.4] model (the mass configuration 
studied by zur Linden 5 ) and also by Kaufmann 7 )) the eigenvalues of the first and second 
excited states form a complex conjugate pair in the mass range 0.25 < s < 2.65. We shall 
analyze here the eigenvalue problem for the p-wave (1 = 1) case by a more accurate numerical 
method than that used in Ref. 8. We have truncated the infinite sum in Eq. (2-3) over V 
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(from 1 to oo) to a finite sum (from 1 to a cutoff value L c ), and replaced the integration over 
z' by the Appoint Gauss-Legendre quadrature formula. By this approximation, the infinite 
system of integral equations (2-3) is reduced to the L c ■ iV-th order real but nonsymmetric 
matrix eigenvalue problem: 

L c N 

9L,i(zj) = A Y K LL',i(zj,z k )w k gv,i{z k ) . 

L'=l k=l 

(L = 1,2,...,L C , j = l,2,---,N) 

The points Z\ , • • • , z^ are the Gauss-Legendre points on the interval [—1, 1], and wi, • • • , wn 
are the corresponding integration weights. 

A preliminary numerical calculation 
made earlier for the p-wave [1.6 <= 1.0 =>- 
0.4] model suggested that all the eigenval- 
ues were real. To judge whether the emer- 
gence of complex eigenvalues is a peculiar 
phenomenon to the s-wave case, it is desir- 
able to investigate the eigenvalues in other 
mass configurations. If we fix A and vary /i, 
we can expect that complex eigenvalues ap- 
pear more easily for small \x than for large [i: 
At /j, — (the Wick-Cutkosky model), eigen- 
values of the first and second excited states 
degenerate (forming multiple poles consist- 
ing of coinciding simple poles in the scatter- 
ing Green's function 10 )) at s = A A 2 (pseu- 
dothreshold). If /i aquires a small nonzero 
value, it is probable that the degeneracy is 
removed and two coinciding simple poles be- 
come a complex conjugate pair for a nonsymmetric real matrix case. For this reason we have 
lowered the value of \x and investigated the eigenvalue problem for the [1.6 <= 0.2 =^ 0.4] 
model. We have calculated first six eigenvalues, and have found that the second and third 
eigenvalues and also the fourth and fifth eigenvalues form complex conjugate pairs for s in 
the range 0.60 < s < 2.25 and 0.48 < s < 2.35, respectively. We have sketched in Fig. 1 
a graph of the first six eigenvalues as a function of s (0 < s < 4) in the complex A-plane. 
At s = 4Z\ 2 , for example, A 2 , 3 = 11.9 =F 0.41i and A 4i5 = 22.6 =F 0.98i The ground state Ai 
is always real since there is no partner coinciding with it at s = AA 2 and /i — 0. Also the 



Re(X) 




s 



Fig. 1. The eigenvalues of unequal- mass scalar- 
scalar model with A = 0.6 and fi = 0.2. N 
(the number of Gauss-Legendre points) is 
47, and L c (the four-dimensional angular 
momentum cutoff) is 10. 
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sixth eigenvalue A6 remains real: The 4-, 5- and 6-th eigenvalues degenerate at s = 4A 2 and 
\x = 0, and the lower two A4 and A5 form a pair, while A6 sticks to a life of proud loneliness. 
We have also confirmed that in the [1.6 -<= 0.2 =>- 0.4] model complex eigenvalues certainly 
appear in the s-wave case. 

§4. BS amplitudes 

We shall investigate the behavior of the BS amplitude from the numerically obtained 
p-wave eigenvectors g^\(z; s), corresponding to the k-th eigenvalue A = A& (k — 1, 2, • • • , 6) 
at the bound-state mass squared s. The radial part (with the factor Yi m (9,(f)) omitted) of 
the approximate BS amplitude <f>k(\p\,P4] s ) is given by 

*(|p|.M«) = |p||/») I - 1 ( r ^) I+3 ^(^) S S(^; S ) ■ (4-i) 

As the first example, we shall consider a rather deep relativistic case s = 0.40, where all 
the eigenvalues are real: Ai = 6.45, A 2 = 14.1, A 3 = 14.7, A 4 = 26.2, A 5 = 26.7, A 6 = 27.7. 
We can see that A2 and A3 are nearly equal, so are the A4, A5 and A6- This can be understood 
as a remnant of the degeneracy between A2 and A3 (and also among A4, A5 and A6 ) at the 
pseudothreshold in the zero-mass exchange case. We have plotted in Figs. 2(a) and 2(b) the 
real and imaginary parts of the BS amplitude of the first excited state {k = 2) on the P4-\p\ 




Fig. 2. The BS amplitude of the p-wave first excited state (corresponding to the real eigenvalue) 
at s=0.40, (a) real part and (b) imaginary part. 

half plane. The complex-valued amplitude is normalized so that at the point where the 
absolute value becomes largest, the real part of the magnitude is equal to 1.00, while its 
imaginary part is zero. We can see that due to the presence of the kinematical factor \p\ 
in Eq. (4T), the BS amplitude vanishes on the p 4 -axis and peak of the amplitude is pushed 
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into the inner |p] region. The real (imaginary) part is an even (odd) function of the variable 
Pa- 

The real and imaginary parts of the BS amplitude corresponding the complex eigenvalue 
A2 = 11.9 — 0.41i at s = 1.44 (pseudothreshold) are depicted in Figs. 3(a) and 3(b), respec- 
tively. Other eigenvalues are Ai = 5.1, A^s = 22.6 =F 0.982 and X e = 23.1. The contour plot 
on the horizontal plane shows that the symmetry (antisymmetry) of the real (imaginary) 
part of the amplitude under the ^-inversion is broken. 




Fig. 3. The BS amplitude of the p-wave first excited state (corresponding to the complex eigen- 
value) at s=1.44, (a) real part and (b) imaginary part. 

As the last example, the BS amplitude of the ground state at the nearly nonrelativistic 
energy s = 3.90 is plotted in Fig. 4. The corresponding eigenvalue is Ai = 1.2. The eigenvalue 



1 n 0.04- 




Fig. 4. The BS amplitude of the p-wave ground state at s=3.90, (a) real part and (b) imaginary 
part. Note that the scale of the vertical axis in Fig. 4(b) is compressed to one tenth compared 
with that in Fig. 3(b). 

approaches zero in the nonrelativistic limit (s j 4), since the ground state have a well- 
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defined limit (normal state). At this energy, the ^-symmetric real part dominates over the 
P4-antisymmetric imaginary part: The maximum of the imaginary part is at most several 
percent of the maximum of the real part. Moreover, if we integrate the amplitude over the P4 
variable, which corresponds to the equal-time nonrelativistic wave function, the contribution 
from the imaginary part vanishes identically. As in the s-wave case, the support of the 
amplitude is concentrated near the origin. 

§5. Summary and discussion 

We have thus confirmed that the emergence of complex eigenvalues (of the coupling 
constant) is a commonplace event in the scalar-scalar ladder model in the unequal-mass case. 
For a fixed mass difference 2 A between two constituent paticles, complex eigenvalues in the p- 
wave case appear in the narrower mass range of the exchanged particle than that in the s-wave 
case. It is expected that complex eigenvalues also appear in higher waves if the exchanged- 
particle mass ji is nonzero and sufficiently small. This is due to the degeneracy of eigenvalues 
of the Wick-Cutkosky model at the pseudothreshold s = A A 2 : A pair of eigenvalues, which 
is degenerate at s = 4Z\ 2 when /i — 0, is attracted into a complex conjugate pair when 
fi becomes nonzero, because of the non-Hermiticity of the matrix equation when A is not 
equal to zero. Based on these considerations, it will be possible to prove analytically (or 
at least semi- analytically) the existence of complex eigenvalues at s = 4Z\ 2 , with the aid of 
perturbation theory around fi = 0. 
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